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Although the effects of interactions in solid state systems still remains a widely 
open subject, some limiting cases such as the three dimensional Fermi liquid or the 
one-dimensional Luttinger liquid are by now well understood when one is dealing 
with interacting electrons in periodic crystalline structures. This problem is much 
more fascinating when periodicity is lacking as it is the case in quasicrystalline 
structures. Here, we discuss the influence of the interactions in quasicrystals and 
show, on a controlled one-dimensional model, that they lead to anomalous trans- 
port properties, intermediate between those of an interacting electron gas in a 
periodic and in a disordered potential. 

1 Introduction 

Our progress in solid state physics have been drastically dependent on our under- 
standing of the physical properties of periodic structures. The knowledge of the 
eigenstates (Bloch waves) in such structures, enabled us to tackle other physical 
effects such as electron-electron interactions and transport properties. Some stud- 
ies have also been made to understand systems totally lacking periodicity such as 
disordered ones. In this case also the knowledge of the diffusive, or exponentially 
localized nature of the wavefunctions has been a good starting point. An intriguing 
situation occurs in the case of quasicrystals H, for which the system is neither peri- 
odic nor disordered. On the experimental side the transport exhibits many unusual 
properties. These metallic alloys are notably characterized by a low electrical con- 
ductivity a which increases when either temperature or disorder increases The 
very low temperature behaviour of a is still an open question, and strongly depends 
on the materials. For example, in AlCuFe u and AlCuRu a finite conductivity 
at zero temperature is expected whereas recent results til seem to confirm a Mott's 
variable range hopping mechanism (cr(T) ~ exp — (Tq/T) 1 ^ 4 ) for i — AlPdRe icosa- 
hedral phase down to 20 mK. The optical conductivity is also unusual since there 
is no Drude Peak for icosahedral quasicrystals Bu. 

From a theoretical point of view, the case of independent electrons in one- 
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dimensional (ID) systems has been deeply investigated for different quasiperiodic 
structures (Harper model, Fibonacci chain, ..)&□, giving rise to singular continuous 
spectra with an infinite number of gaps. Moreover, the corresponding eigenstates 
are neither extended not. .localized but critical, and are known to be responsible 
of anomalous diffusion liiM. For higher dimensional systems (Penrose tiling, Oc- 
tagonal tiling, icosahedral structure...), similar studies had also displayed coraclcx 
and intricated spectra, with analogous characteristics of the electronic states Ii3£j. 
Finally, studies of transport properties^SESOlla has also displayed anomalous 
aviour such as a power-law decreasing of the conductance with the system size 



hglt 

Given the difficulty of the full quasiperiodic problem, and the relatively poor 
knowledge of the single particle wavefunctions, adding interactions is a very diffi- 
^Iven in ID incommensurate structures, few results have been obtained 
In order to treat interactions in quasicrystals we thus follow a different 
route, used with success for disordered systems. We first solve the periodic system 
in presence of interactions. This is relatively easy, either in the onc-dimcnsional 
case which is considered in this paper, for which exact solutions exists, or even in 
higher dimensions through approximate (Fermi liquid) solutions. Then, we study 
the effect of a weak perturbative quasiperiodic potential on this solution via a renor- 
malization group approach. The main idea consists in analyzing the dependence of 
the low-energy physics with respect to the Fourier spectrum of the potential. The 
influence of the renormalization group procedure is discussed in terms of an effec- 
tive elimination of the various potential harmonics around the Fermi leveL leading 
to unusual properties in the case of a self-similar quasiperiodic structureES. 

2 The model 

Let us consider a system of interacting spinless fermions on a lattice embedded in 
a general on-site (diagonal) potential W, described by the following Hamiltonian: 

H = -t^c\c J +V^n l n l+l +^2w i n l (1) 
<i,i) 1 1 

= H (t,V) + H w (2) 

where cj (resp. Cj) denotes the creation (resp. annihilation) fermion operator, 
rii = cj Ci represents the fermion density on site i, and (. . .} stands for nearest 
neighbors pairs. ._. 
In the continuum limit, Ho(t : V) writeso: 

Ho(t, V)^^jdx [(uK)(nU) 2 + (c^) 2 ] (3) 

where is a boson field related to the long wave length part of the fermionic density 
by p(x) = ~V 4>(x) / it 1 and II is its canonically conjugate field. This approximation 
provides the correct description of the low energy physics. All the interactions are 
absorbed in the two constants u and K. u is the renormalizcd Fermi velocity (in 
the non interacting case: u = vp = 2iasin(fcira)). K is the Luttinger parameter 
controlling the decay of various correlation functions: K=l in the non-interacting 
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case, K > 1 for attractive interactions (V < 0) and K < 1 for repulsive interactions 
(V > 0). Note that, for many models..tbese parameters can be explicitely computed 
in terms of the microscopic details Eal2j. 

The perturbative part of H , can also be expressed in terms of these boson fields: 

H w = — [ dx W(x) cos [2k F x - 2(/>(x)] (4) 
na J 

where a is a short distance cut-off of the order of the lattice constant a. 
3 Renormalization group analysis 

To determine the long distance physics, we use a standard, perturbative renormal- 
ization group approach. Full details can be found in E£rE3 and we focuss here on 
the main equations for the potential and interaction parameter: 

f - 4<*> 



with 

G(I)=^^^J[( ? I*)«(I)]. (6) 

E=±l q 

y q = aW(q)/u is the dimensionless Fourier components of W and a(l) = a(0) e l 
is the renormalized short distance cutoff. In Eq. (^J), J is an ultraviolet regulator 
whose precise form depends on the procedure used to eliminate the high energy 
degrees of freedom. Typically one has: 

J{x) ~ 1 if x < 1 (7) 
= otherwise (8) 

Physically, this renormalization is equivalent to an investigation of the low- 
energy properties in a window around 2k F in the reciprocal space whose width is 
proportionnal to e~ l . Thus, the full Fourier landscape of the potential determine 
the scaling of the parameters u and and one has to distinguish several situations. 

If the Fourier spectrum of W has only a single peak or few well separated peaks, 
then the system is analogous to a periodic system. If the Fermi level is not on one of 
the peaks then at some lengthscale the regulator J gives always zero (see figure 0) . 
The potential W is irrelevant and one recovers a Luttinger liquid metallic phase. 
On the contrary, if the Fermi level is right on one of the peaks (e.g. 2k F = qo), all 
physical properties are dominated by this single harmonic (the regulator cuts all 
the others) and the system feel a periodic-potential W(x) = Xcos(qox). The case 
of the single harmonic is well known HEBEj: the periodic potential is relevant for 
K < 2 and opens a gap 

A ~ yf* (9) 

For K > 2 the periodic potential is irrelevant and system remains a Luttinger liquid 
with gapless excitations. Our equations clearly show that novel physics occurs when 
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(b) 



Ef q 




Figure 1. Fourier spectrum of the potential W for the periodic (a), disordered (b) or quasipcr- 
diodic (c) cases. All the peaks within a window of width ~ e~ 1 (see text) around the Fermi energy 
Ep (shown as a box) control the physics of the problem. This is crucial for the disordered (b) 
and quasipcriodic (c) systems for which it is impossible to really isolate a single peak, contrarily 
to the case of periodic systems (a) 
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the spectrum of the potential W is dense. This is the case both for disordered 
systems and for the quasiperiodic ones (see figure [l]). Indeed in that case it is never 
possible to really isolate a single peak O. The physics is thus not controlled by 
the presence or not of a peak at the Fermi level but by the forest of peaks in the 
(shrinking) window of width e~ l . In this respect, the properties of a quasicrystal 
are much closer (with proper differences discussed below) to the one of a disordered 
system than the one of a "periodic" (i.e. with a limited number of harmonics) one. 

The disordered case can easily be obtained from our equations. Let us assume 
that W is uncorrelated with averages given by 

W*(q)W(q') = A 2 V (10) 

In the limit of weak disorder, (g) can be integrated neglecting the renormalization 
of K: y(l) = y(0)e {2 - K ^ 1 . Then, (g) simply becomes: 

^ = -if 2 C e < 3 - 2if > ! (11) 

where C is a constant. Eq. (|Tl| ) provides a critical value K c = 3/2, separating 
an Anderson-like insulating phase (K < K c ) from a metallic state (K > K c ) 
Note that contrary to the periodic case, this metal-insulator transition occurs for 
any value of kp. Indeed, as shown figure 0, the Fourier spectrum is uniform (after 
averaging) so that all the positions of the Fermi level are equivalent. 



4 The quasiperiodic case 

The quasiperiodic case provides an interesting intermediate situation. Let us con- 
sider for simplicity a Fibonacci potential although our results can be easily extended 
to more general situations. 

The quasiperiodicity is provided by the Wi's that take two discrete values Wa = 
+A/2 or Wb = —A/2 given by the spatial modulation of the Fibonacci chain. In 
fact, we consider a periodic approximant of this structure with Fj sites per unit cell 
that can be obtained by j iterations of the substitution rules: A — » AB, B — > A, 
where Fj is the j element of the Fibonacci sequence defined by: 

Fi = F 2 = 1 

F j+1 =F j +F j _ 1 (12) 

We denote p — i^-_2,s = Fj—x,n — Fj and n! — s (resp. n! = p) if j is even 
(resp. odd). In the quasiperiodic limit (j — > oo), the ratio s/p converges toward 
the golden mean r = ■ Tha.Fourier transform of W can be straightforwardly 
using the conumbering scheme E3: 



2irm 



W [9= = . (nmn \ (13) 

for to = 1 to n — 1 (a is the lattice spacing). A global shift of the Wi allows us to 
deal with a zero-averaged potential so that we can set W(0) = 0. Moreover, the 
underlying substitution rule provides a self-similar structure that can be readily 
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seen m Fig. In the non interacting case, each Fourier component of the potential 
opens a gap whose width is giyen, at first order in perturbation, by the amplitude of 
the corresponding componentEj. In the quasiperiodic limit, the Fourier spectrum of 
the potential become dense so that the spectral measure goes to zero. Therefore, for 
any filling factor, the Fermi velocity vanishes and the system is an insulator (at zero 
temperature). In contrast, given the complexity of the Fourier spectrum, one can 
expect the long distance physics of the interacting system to depend on the position 
of the Fermi level. Indeed, the flow equations (|J) show that the behaviour of G 
is strongly sensitive to the position of the window in the spectrum. For a given 
maximum renormalization length scale l m ax , corresponding to accessible physical 
range, two different cases must be distinguished. 

First, if the Fermi momentum 2fcp is close to a main peak of the Fourier spec- 
trum, i. e. (q ± 2k f ) _1 > lmax then, at long distance, i.e. up to I ~ lmax, the flow 
of K is controlled by this harmonic (see Fig. [j] (upper curve)). A metal- insulator 
transition can occurs with K c = 2. This case is very similar to the simple periodic 
one. 

A more interesting behaviour is encountered when the Fermi level is far from a 
dominant harmonic of the quasiperiodic potential, for example at half-filling [2kp = 
ir /a). Indeed, the low energy properties up to l ma x are no more dominated by the 
ultimate presence of a gap or not but by the precise dependence of G with the scale. 
As can be seen from figure G has an exponential scaling G ~ e^ 2 ~ 2K H. Such a 




Figure 2. Scaling of the function G. The presence of many peaks gives a scaling dimension 
different both from disordered or periodic systems. 
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behaviour has several consequences. First it provides a critical value of K = 1 
for the metal-insulator transition. Let us point out that this value is the smallest 
critical value that one can expect for such a transition. Indeed it is clear that 
repulsive interactions (K < 1) enforce the insulating character of the system. To 
know whether this analysis is asymptotically correct or not, one would clearly need 
an analytical calculation of G, a rather complicated task. Note however that our 
predictions seem to be confirmed by numerical work based on the Density Matrix 
Renormalization Group method on the XXZ chain with quasiperiodic exchange Ea. 
This work provides evidence for different behaviours for the scaling of the gap with 
the size for various K. In particular, a transition is observed for K = 1 at half-filling 
even in the strong coupling regime. An interesting, and yet unanswered question is 
how the spectrum, or more generally the density of states is modified when K varies. 
It is very natural to assume that the small gaps close up when K increases, given 
that the system is always gapless for K — 2. The spectrum would thus evolve from 
a set of zero measure for K — 1 to a set of finite measure for K > 1 as indicated in 
Figure |3|. It would be interesting to check this scenario on the XXZ chains under a 
magnetic field. So, the quasicrystal differs from a periodic one, for which the gap 

— • 1 1 1 ■ — K=l 



1 <K<2 



K=2 

Figure 3. A possible evolution of the density of states as a function of the Luttinger liquid 
parameter K. For K = 1 the spectrum is a set of measure zero and gaps are present almost 
everywhere. For 1 < K < 2 the small gaps close as K increases. For K = 2 even the largest gaps 
are closed since even a periodic potential is irrelevant above K = 2. 

only acts for a given position of the Fermi level with K c — 2, and a disordered 
system for which the potential is relevant regardless of the position of the Fermi 
level, but below a constant critical value K c = 3/2. This important modification 
of K c is reminiscent of a correlated disorder with long range correlations in space 

for which the averaged disorder potential W*(q)W(q') = S qq >A(q) is not constant. 

The other important consequences of such a behaviour are of course exhibited 
by the transport properties. A straightforward calculation would give a power law 
for the frequency dependence of the conductivity cr(w) ~ (l/ui) 5 ~ 2K ~ t * with fx = 2 
(at variance with \i = for periodic and /j, = 1 for disordered). This anomalous 
power law depending both on the interactions and the spectrum properties would 
also appear in a size dependence of the Landauer conductance R{L) ~ L 4_2A_Al . 



quasi '06: submitted to World Scientific on February 1, 2008 7 



The temperature dependence of the conductivity raises an interesting question, 
similar to the one occuring for periodic systems. Cutting the flow at a the natural 
lengthscale / = log(vp/T) would lead to a power law dependence of the conductivity 
cr(T) ~ j'3-2-K"-At j anc j t nus to a finite conductivity. In order to obtain such a 
behaviour it is probably necessary to introduce dephasing processes (like a coupling 
to a thermal bath). If such a temperature dependence of the conductivity occurs 
it would be again a demonstration of the close connection between quasiperiodic 
systems and disordered one (for which it is quite natural to expect dissipation at 
finite temperature). 

Let us note that is it possible to investigate the effects of disorder on the 
quasiperiodic potential using (||). Disorder in d = 1 leads in itself to localization, 
which is to be avoided if one wants to make connection with higher dimensions. 
One way to achieve this is simply— to introduce forward scattering on impurities 
which cannot lead to localization E3. Introducing such a disorder cuts the flow in 
at a characteristic lengthscale £f inversely proportional to the disorder. Since it 
prevents the (quasi-)periodic potentiallo grow it makes thus the system more con- 
ducting (more details can be found inEJ). One thus recovers on this simple model 
the physical effects leading to the inverse Mathiessen law experimentally observed. 

Exploring these issues, as well as investigating the consequences in higher di- 
mension is a challenging problem left for the future. 



References 

1. D. Shechtman, I. Blcch, D. Gratias, and J. W. Cahn. Phys. Rev. Lett., 
53:1951, 1984. 

2. D. Mayou, C. Berger, F. Cyrot-Lackmann, T. Klein, and P.Lanco. Phys. Rev. 
Lett, 70:3915, 1993. 

3. T. Klein, C. Berger, D. Mayou, and F. Cyrot-Lackmann. Phys. Rev. Lett., 
66:2907, 1991. 

4. B. D. Biggs, S. J. Poon, and N. Munirathan. Phys. Rev. Lett, 65:2700, 1990. 

5. J. Delahaye, J. P. Brison, and C. Berger. Phys. Rev. Lett., 81:4204, 1998. 

6. C. Homes, T. Timusk, X. Wu, Z. Altounian, A. Sahnoune, and J. O. Strom- 
Olsen. Phys. Rev. Lett, 67:2694, 1991. 

7. S. Burkov, T. Timusk, and N. W. Aschroft. J. Phys. C, 4:9447, 1992. 

8. M. Kohmoto, L. P. Kadanoff, and C. Tang. Phys. Rev. Lett, 50:1870, 1983. 

9. S. Ostlund, R. Pandit, D. Rand, H. J. Schellnhuber, and E. D. Siggia. Phys. 
Rev. Lett, 50:1873, 1983. 

10. J. X. Zhong and R. Mosseri. J. Phys. C, 7:8383, 1995. 

11. F. Piechon. Phys. Rev. Lett., 76:4375, 1996. 

12. B. Passaro, C. Sire, and V. G. Benza. Phys. Rev. B, 46:13751, 1992. 

13. S. Yamamoto and T. Fujiwara. Phys. Rev. B, 51:8841, 1995. 



14. D. Mayou, 1999, =ond-mat 9912113 



15. C. Sire. Lectures on Quasicrystals. Editions de Physique, Les Ulis France, 
1994. 

16. S. Roche, D. Mayou, and G. Trambly de Laissardiere. J. Math. Phys., 38:1794, 
1997. 



quasi '06: submitted to World Scientific on February 1, 2008 



17. S. Roche and T. Fujiwara. Phys. Rev. B, 58:11338, 1998. 

18. H. Schulz-Baldes and J. Bellissard. J. Stat. Phys., 91:991, 1998. 

19. H. Tsunctsugu and K. Ueda. Phys. Rev. B, 43:8892, 1991. 

20. Masaki Goda and Haruhiko Kubo. J. Phys. Soc. Jpn., 58:2109, 1989. 

21. H. Hiramoto. J. Phys. Soc. Jpn., 59:811, 1990. 

22. J. C. Chaves and I. I. Satija. Phys Rev. B, 55:14076 1997. 

23. J. C. Chaves and I. I. Sat ija, 1998. |cond-mat 9803104 

24. Vieri Mastropietro, 1998 |cond-mat 981012Sf 



25. D. Sen and S. Lai, 1998. cond-mat 9811330 



26. J. Vidal, D. Mouhanna, and T. Giamarchi. Phys. Rev. Lett., 83:3908, 1999. 

27. H. J. Schulz. Elsevier, Amsterdam, 1995. 

28. F. D. M. Haldane. J. Phys. C, 14:2585, 1981. 

29. F. D. M. Haldane. Phys. Rev. Lett, 45:1358, 1980. 

30. J. Vidal et al., to be published. 

31. A. Luther and V. J. Emery. Phys. Rev. Lett, 33:589, 1974. 

32. T. Giamarchi. Phys. Rev. B, 44:2905, 1991. 

33. E. Kolomeisky. Phys. Rev. B, 47:6193, 1993. 

34. This point was missed in previous studies p£ a quasiperiodic potential, that 
considered a Harper like potential cos(Qx) Indeed, as was already known 
for the noninteracting case, such a potential is perturbatively very similar 
to a simple periodic one, and quasiperiodic effects only occur markedly for 
large enough potential, making it unsuitable for an RG study. The Fibonacci 
potential studied here does not suffer from such problems and exhibit the 
quasiperiodic effects even at small coupling. 

35. T. Giamarchi and H. J. Schulz. Phys. Rev. B, 37:325, 1988. 

36. W. Apel. J. Phys. C, 15:1973, 1982. 

37. R. Mosscri. In J. Yacaman et al., editor, Proceedings of the 3rd International 
Meeting on Quasicrystals, page 129, Singapore, 1990. World Scientific. 

38. C. Sire and R. Mosseri. J. Phys. (Pans), 50:3447, 1989. 

39. K. Hida, 1999. cond-mat 9911424. 



quasi '06: submitted to World Scientific on February 1, 2008 



